Exploring spatio-temporal changes in coastal recreational fisheries and potential links to extreme weather events

Extreme weather events across coastal environments are expected to increase in frequency under predicted climate change scenarios. These events can impact coastal recreational fisheries and their supporting ecosystems by influencing the productivity of fish stocks or altering behaviours and decision-making among fishers. Using off-site telephone/diary survey data on estuarine and oceanic recreational fishing activity in eastern Australia, we analyse interannual and geographic variability in bream (Acanthopagrus spp) and snapper (Chrysophrys auratus) catch, total effort and total catch per unit effort (CPUE) through a period (2013/2014, 2017/2018 and 2019/2020) that encompassed severe drought, bushfires and flooding. Interacting spatial and temporal differences were detected for bream and may reflect spatial variation in the intensity and extent of some of the extreme weather events. The catch of snapper did not change temporally, providing little evidence that this species’ catch may be influenced by the extreme weather events. Independent bioregional and temporal effects on effort were detected, while CPUE only showed significant bioregional differences. Although adverse conditions created by the extreme weather events may have dissuaded fisher participation and impacted effort, we propose that the observed temporal patterns in effort reflect the early influence of socio-economic changes brought on by the COVID-19 pandemic on coastal recreational fishing, over and above the impacts of extreme weather events. This study demonstrates how interrelated ecological, social and economic factors can shape coastal recreational fisheries and facilitates development of management strategies to address future threats to the sector.


Introduction
The frequency and intensity of extreme weather events such as hurricanes, heat waves, drought, floods and wildfires and is predicted to increase under projections of climate change [1,2].Although direct attribution of extreme weather events to anthropogenic climate change is a topic of ongoing scientific debate [3,4], the trend of increased extreme events has recently become apparent in many parts of the world [5][6][7].These events can have significant and measurable impacts on coastal recreational fisheries, which shape the ecology and play vital socioeconomic roles in the marine resource landscape of many high-income and transitional nations [8,9].Like many other factors, extreme weather events can influence the productivity of fish stocks or alter behaviours and decision-making among fishers and thereby impact recreational fishery metrics such as catch and catch per unit effort (CPUE) [10,11].
Coastal recreational fisheries can be directly or indirectly affected by extreme weather events.Direct effects include preclusion of fishery access through bad weather or damage to fishing infrastructure such as shorelines, marinas, boats and fishing gear; as well as loss of fish biomass due to displacement or mortality.For instance, an extreme cold event in southern Florida during 2010 resulted in widespread fish-kills for the common snook, Centropomus undecimalis [12].This species supports a substantial coastal recreational fishery that experienced subsequent reductions in catch per unit effort within western estuaries [12].Indirect effects of extreme weather events on coastal recreational fisheries can include the loss of coastal habitats such as marshes, oyster reefs and mangroves which support the productivity of exploited fish populations.Duke, Kovacs [13] described how extreme high temperatures and low precipitation contributed to extensive dieback of mangrove tidal vegetation in Australia's Gulf of Carpentaria and inferred a reduction in harvest for recreational fisheries targeting marine finfish (e.g., king salmon (Polydactylus macrochir)) that depend on mangrove habitats.
Monitoring recreational fisheries and integrating dynamic management actions in response to extreme weather events can provide social and economic benefits to resource users, support ecosystem benefits and promote fisheries sustainability [14,15].Recreational fishing surveys are useful tools for monitoring temporal variation in fishery indices, particularly as the balance of effort in fisheries shifts among sectors and recreational catches of many species exceed that of other sectors.However, the large, diffuse and heterogenous nature of many recreational fisheries can make it difficult to collect reliable and representative on-site data in a cost-effective way [16].Offsite telephone/diary surveys based on randomized sampling overcome many of these challenges by enabling broad yet efficient fishery coverage due to larger sample sizes and/or higher response rates [17][18][19].Such surveys are consequently used to routinely monitor and assess expansive recreational fisheries in nations such as the USA, Germany and Australia [20][21][22][23].
In the Australian state of New South Wales (NSW), off-site telephone/diary surveys have been implemented over the last decade to track spatio-temporal patterns of recreational fishing activity [24,25].NSW is home to the nation's greatest number of recreational fishers, with participation rates of approximately 11.9% [21,26].The coastal saltwater fishery consists of estuarine and oceanic components across 2,000 km of coastline and is the largest in the state [26,27].The recreational fisheries of NSW have been monitored using a standardised annual survey design applied only in 2013/14; 2017/18 and 2019/20.The last two surveys coincided with periods during which three extreme weather events occurred in NSW, potentially impacting coastal aquatic environments [28,29]: (i) the most severe drought recorded in the European history of Australia [30]; (ii) the "Black Summer" bushfires which burnt through more land than any fires in the past 25 years [5,6]; and, (iii) heavy rainfall and subsequent flooding that led to some of the highest river levels since 1992 [7].The synchronicity of these extreme weather events with the two later surveys provided a unique opportunity to examine changes in annual recreational fishery metrics in relation to the extreme weather events.
This study therefore reports on temporal patterns in coastal recreational fishing activity through 2013/14, 2017/18 and 2019/20 to understand shifts in the fishery that coincided with the extreme weather events.Focusing on the line-based recreational fishery (also known as angling) of NSW, which represents the dominant (97% of participation) form of fishing [24,31], this study specifically compares estuarine and oceanic fishing effort, CPUE and the species-specific total catch of Acanthopagrus spp (species complex of A. butcheri [black bream] and A. australis [yellowfin bream], hereafter referred to as bream) and Chrysophrys auratus (hereafter referred to as snapper) across the study period.These species are consistently among the five most commonly caught species in both the estuarine and oceanic waters of NSW [24,31] and are also of significant commercial importance within the state, contributing ~$2M per annum in production value [32][33][34].Furthermore, both species are predominantly captured using line-based fishing in NSW [24,31].Since the extent of the extreme weather events in question was not spatially homogenous and recreational activity in NSW has been shown to vary spatially [25,35], this study also examines whether coarse spatial interactions occurred with temporal patterns of recreational fishing.This research explores how interrelated ecological, social and economic drivers can potentially shape coastal marine recreational fisheries in a rapidly changing world.

Materials and methods
This study, involving human participants, was approved by the New South Wales Department of Primary Industries Human Ethics process (INT20/76587).

Study area
This study considered estuarine and oceanic recreational fishing data collected across the NSW coastline (also known as the NSW Marine Estate; NSW Marine Estate Management Authority [36]).This coastline is divided into 3 geographic management bioregions by the state's Marine Estate Management Authority: the North Coast (which extends from the Queensland border (28.16˚S; 153.55˚E) down to Stockton Beach (32.81˚S; 151.96˚E)); the Central Coast (which extends from Stockton Beach to Shellharbour (34.58˚S; 150.87˚E)); and the South Coast (which extends from Shellharbour to the Victorian boarder (37.51˚S; 149.98E )) [36,37] (Fig 1).For this study, estuarine waters were defined as saltwater estuaries, bays and inlets while oceanic waters were those that spanned from the coastline to the 200 nm exclusive economic zone (Fig 1).In this context, coastal waters encompass both estuarine and oceanic environments.The 3 geographic management areas were affected by severe drought, bushfires and flooding to varying degrees during the 2013/14, 2017/18 and 2019/20 periods.Table 1 presents a qualitative comparison of the relative intensity of each of these extreme weather events in each geographic area in each year.Although the COVID-19 pandemic was not an extreme weather event of initial interest in this study, it was an unprecedented socioeconomic shock that took place during the 2019/20 and, as discussed later, had the potential to influence the fishery metrics of interest as was reported in other parts of the world [38,39].Its relative intensity is therefore also depicted in Table 1.

Data collection
Annual data on the recreational fishery was collected for fishing activity that occurred during the 2013/14, 2017/18 and 2019/20 periods.For each period, data was collected using an off-site telephone-diary survey developed to provide cost-effective state-wide fishery information over a large spatial scale [21,26,49].The approach involved a two-phase design, with an initial Screening Phase and the second an intensive Diary Phase.The Screening Phase occurred from March to May 2013, September to October 2017 and August to October 2019.It was a structured telephone interview of a randomly selected sample of households from the NSW  Recreational Fishing Licence (RFL) database of 1-3 year licence holders [31].During the screening phase, profiling information was collected (including age, gender and fishing avidity) and respondents answered questions regarding their intention to fish in the ensuing 12 months to determine their eligibility for the Diary Phase.A longitudinal panel survey made up the study's Diary Phase and monitored the fishing activity of all residents (aged � 5 years) within recruited households between June 2013 to May 2014, between October 2017 to September 2018 and between November 2019 and October 2020 [24,31].Fishing information collected during the Diary Phase included date, location, start and finishing times, fishing methods, fishing platform (boat-or shore-based), waterbody type (freshwater, estuarine or oceanic), target species and numbers of kept and released animals by species/species group [31].This information was recorded by trained interviewers during regular telephone calls to diarists (which were made on at least a monthly basis to reduce recall bias) [31].The Screening Phase and Diary Phase are described in greater detail in Ochwada-Doyle, Miles [50] and Murphy, Ochwada-Doyle [31].

Data analysis
The statistical programming software R [51] was used for all analyses and graphing described hereafter.To enable estimation of total and mean population catch, effort, and CPUE, analysis of each year's raw survey data was based on a stratified random survey design.The RFL holders (n) represented the primary sampling unit within a randomly selected household, and the fishing activity of all resident fishers within the household was surveyed (representing a form of single stage cluster sampling).Adjustment and calibration were integrated to expand information from each sampled household to population estimates.A response propensity model was fitted to partially-responding (assumed to be representative of non-responders) and fullyresponding RFL holders who provided information on the fishing activity of household residents in the last 12 months.This model (which predicted a fisher's unique combination of characteristics that are expected to affect the relative likelihood of obtaining a complete survey response from that fisher [52,53]) was then used to adjust sample data for response bias [26,53].Calibration (weighting adjustment to make the estimated and known population totals consistent with each other) also implemented using benchmark data (known population total of RFL holders) containing the total number of licenced holders in each of the state's ten residential strata [26,53,54].In order to account for unexpected "drop-ins" into the fishery, an additional adjustment was applied based on a "non-intending call-back" survey.This accounted for the avidity index (which was measured for household residents aged �5 yrs.as a function of estimated fishing frequency in the previous 12 months [53]) reported for "dropins" and residential stratum [26,53].The Survey [54] package was used to expand and calibrate the raw data and provide total population estimates of species-specific catch and total combined effort following the approach outlined in Lyle et al. [53] (see Lyle, Wotherspoon and Stark [53] and Lumley [54] for a full description and equations).Mean population estimates of species-specific catch and total effort were also generated from the expanded raw data on species-specific catch and total combined effort.The raw data on total catch across all species and total combined effort were used to calculate total CPUE as a function of numbers caught/fisher day for each fisher.They were then expanded and used to estimate mean CPUE across RFL households.
For catch, effort and CPUE data from each survey period, separate analyses were conducted for data from oceanic and estuarine environments.Using the glmer.nb()function in the lme4 package [55], generalized linear mixed-effects models (GLMMs) assuming a negative binomial distribution [56,57] were used to examine whether the mean of expanded estimates of species-specific line-based catch (across retained and released catch) for each of bream and snapper (in terms of numbers) changed among the three survey years (2013/14, 2017/18 and 2019/20) and bioregions (North Coast, Central Coast and South Coast).The same analyses were used to examine whether mean expanded line-based effort (in terms of fisher days) changed among the three survey years and bioregions.It's important to note that robust measures of effort could only be estimated across all species captured within a given fishing method due to anglers commonly targeting more than one species during a single fishing event or not specifying a target species in the first place.Similarly, mean expanded CPUE could only be estimated and analysed across all species within a given fishing method.There have been recent attempts to estimate species-specific CPUE from this data which have focused on fishing events where a fisher has identified the species of interest as a target species [58,59].Although this approach can be useful for generating broad performance indicators for a species' fishery, it was not applied in the current study because it relies on the assumption that fishers report their targeting behaviour accurately after their fishing trip (i.e. the reporting is likely to be biased by fishers posteriorly identifying a species as a target species simply because they happened to catch that species in large numbers) and leads to omission of data from trips where there were no specific targets (27% of fishing trips averaged across survey years).In each of the estuarine and oceanic data sets, the line-based CPUE data (in terms of numbers caught across all species per day) had many zeros (>50% of records).The CPUE data were consequently modelled using zero-inflated Poisson GLMMs executed using the mixed_ model() function in the GLMMadaptive package [57,60,61].The use of GLMMs enabled us to account for any non-independence among primary sampling units through inclusion of random-effects terms for individual persons and households [62,63].
The full GLMMs for catch and effort considered all independent parameters and interaction terms and took the general form of: where μ i equals the expected value E(Y i ) and Y i represents the observed catch or effort of a sampled angler after it has been expanded to represent the catch of all anglers in the known population represented by that sampled angler; β 0 is the vertical intercept; β 1 and β 2 are the partial regression coefficients of the regressors for the categorical independent parameters year (x 1,ij ) and bioregion (x 2,ij ) respectively; β 3 is the partial regression coefficient for their two-way interaction; a 1,i and a 2,j respectively represent the random variables person and household; and, ε ij represents the error term [56,64,65].The full zero-inflated GLMMs for CPUE considered all independent parameters and interactions and took on the general form of: In Eq 2, μ i (1-π) equals the expected value E(Y i ) and Y i represents the observed non-zero CPUE value of a sampled angler after it has been expanded to represent the CPUE of all anglers in the known population represented by that angler.β 0-3 , x 1,ij , x 2,ij , a 1,i , a 2,j and ε ij in the first part of Eq 2 (the count component) denote the same attributes as described for Eq 1.The second part of Eq 2 (the zero-inflated component) represents the probability of excess zeros in the data.Here π is the probability of observing a zero; logit() is the inverse logit function that maps the probabilities to the log-odds scale; γ 0 is the vertical intercept; γ 1 and γ 2 are the partial regression coefficients of the regressors for the independent parameters year (x 1,ij ) and bioregion (x 2,ij ) respectively; γ 3 is the partial regression coefficient for their two-way interaction;; v 1,i and v 2,j respectively represent the random effects of person and household on the probability of excess zeros [57,[65][66][67].
Eq 1 and the first part of Eq 2 were initially fitted assuming both a Poisson and a negative binomial distribution family.Akaike Information Criteria (AICs) values were then compared among these models to select the model that had the most appropriate distribution family, whereby the model with the lowest AIC value was deemed most appropriate [64].Based on AIC values, the parameters in the original full model from the most suitable distribution family were iteratively reduced towards parsimony through consecutive exclusion of a parameter [62,68,69].If exclusion of a particular parameter led to a higher AIC value, this meant that the model's parsimony was not improved by that parameter's exclusion and the parameter was therefore retained in the model [64,70].In doing so, the analyses retained the subset of predictor parameters that were most important in explaining variation in catch, effort or CPUE [64].Partial tests have been shown to be an insufficient measure of the appropriateness of alternatives to a model [see 68,71,72], so were only used to assess the influence of an independent parameter or interaction term once the final model had been chosen.For each final parsimonious GLMM, a partial Wald test (α = 0.01) was used to examine the null hypothesis that β i = 0 [55,64,73] for a retained parameter or the interaction term.To enable us to determine the influence of each predictor variable on both the count and zero-inflated components of the final parsimonious zero-inflated GLMMs, partial ɀ-tests (α = 0.01) were used to examine the null hypotheses that β i = 0 or γ i = 0 [55,64,73] for a retained parameter or the interaction term.Where no significant interactions were present and the influence of region was shown to be independently significant, Tukey's post-hoc tests (α = 0.01) were then used to compare the different levels of region in a pair-wise manner with the interaction term from the original model being used as the error term [64].The same procedure was used where no significant interactions were present, and the influence of year was shown to be independently significant.Wherever a parsimonious GLMMs showed that year and bioregion had a significant interaction, a separate secondary series of GLMMs were applied to examine the influence of year at each level of region separately using Wald tests or ɀ-tests (α = 0.01) [55,64,74].This examination was followed by Tukey's post-hoc tests (α = 0.01) whenever year was found to be significant at a particular level of region.

Results
The sample sizes (n = number of RFL holders) used in the survey expansion procedure to estimate the catch, effort or CPUE for a particular year within each bioregion are shown in Table 2.These values also represent the sample size for each year within a region in the GLMMs.

Snapper catch
The most parsimonious GLMM for snapper caught in estuarine waters retained both year and bioregion but excluded their interaction term.Bioregion was the only significant parameter (Table 3) and Tukey's post-hoc tests showed that the total catch of snapper was significantly greater in the Central Coast compared to both the South Coast and North Coast (Table 4; Fig 2A).The most parsimonious model on oceanic snapper catch also retained both independent parameters and excluded their interaction term.The only significant parameter was bioregion (Table 3) and Tukey's post-hoc tests showed that the total oceanic catch of snapper was significantly lower in the South Coast compared to the North Coast and Central Coast (Table 4; Fig 2B).

Total effort
The most parsimonious GLMM for total line-based effort in estuarine waters retained both year and bioregion but excluded their interaction term.Both bioregion and year influenced effort significantly (Table 5).Tukey's post-hoc tests showed that effort, pooled across years, For oceanic data on total effort, the most parsimonious model also retained both independent parameters and excluded their interaction term.Both year and bioregion were significant (Table 5).Tukey's post-hoc tests showed that, spatially, the Central Coast had greater effort than the North Coast and, temporally, 2019/20 had greater effort than 2017/18 (Table 4; Fig 3B).The models, which assumed a negative binomial distribution, applied the Wald test (α = 0.01) to examine the null hypothesis that a β i = 0, where β i s were the regression coefficients for region or year.The total sample sizes (n = number of long-term Recreational Fishing Licence holders) of each model are also shown.

Total CPUE
The most parsimonious zero-inflated GLMM for total line-based CPUE in estuarine waters retained both year and bioregion but excluded their interaction term.Only bioregion had a significant influence on CPUE for both the count and zero-inflated components of the model (Table 6).Tukey's post-hoc tests showed that CPUE, pooled across years, was significantly greater in the Central Coast and the South Coast compared to the North Coast (Table 4; Fig 4A).
The most parsimonious zero-inflated GLMM for total line-based CPUE in oceanic waters excluded year and the interaction term and retained bioregion, which had a significant influence on CPUE for both the count and zero-inflated components of the model (Table 6).Tukey's post-hoc tests showed that CPUE, pooled across years, was significantly greater in the South Coast followed by the North Coast followed by the Central Coast (Table 4;

Discussion
Extreme whether events resulting from increased sea and atmospheric temperatures as well as changes in precipitation can impact coastal recreational fisheries and their supporting ecosystems by influencing the productivity of fish stocks or altering behaviours and decision-making among fishers [10,11].This research monitored interannual variability in bream and snapper catch, total effort and total CPUE across an entire jurisdiction's estuarine and oceanic waters through a period that encompassed severe drought, bushfires and flooding.Each of these extreme weather events impacted different geographic regions and waterbodies variably, with our results suggesting some interactive effects and spatial heterogeneity in terms of fishery metrics.
Discussed below are a selection of the spatio-temporal differences in fishery metrics observed in relation to potential explanatory socio-ecological factors (see Table 7 for summary).While some of the changes in fishery metrics may be linked to the extreme weather events; others are more plausibly explained by considering other factors such as the changes brought about by the COVID-19 pandemic.We stress that there are an array of inter-related mechanisms that could contribute to the fishery patterns observed and therefore pose our explanations as potential links only.To enable less equivocal inferences about the exact causes of the observed spatio-temporal variation, we recommend future experimental work that includes unaffected control sites (for examining localised extreme whether events such floods) that could enable application of a before-after-control-impact (BACI) analytical design [75].A longer time-series (>20 years) of survey data including explicit questions on how often extreme weather events affected a person's choice to fish, angling behaviour and perceived success may also provide robust insights into causes of variation.
The survey showed that for bream caught in estuaries in the North Coast, catches were highest in 2019/20 followed by 2013/14 and lowest in 2017/18 (Table 7; point i).However, the To enable examination of the influence of each predictor variable on both the count and zero-inflated components of each model, partial ɀ-tests (α = 0.01) were used to examine the null hypotheses that β i = 0 or γ i = 0 for a retained parameter or interaction term (where β i and γ i respectively represent the partial regression coefficients for a variable's regressor in the count and zero-inflated components of the model).The total sample sizes (n = number of long-term Recreational Fishing Licence holders) of each model are also shown.number of bream caught in estuaries in the South Coast was significantly lower during 2019/ 20 (Table 7; point i).One explanation for these spatio-temporal differences in the estuarine catch of bream may be linked to disparities in the intensity and extent of the bushfires.The bushfires of the austral spring and summer 2019/20 affected NSW's coastal regions variably, with the South Coast showing the greatest quantity of canopy affected by fire from the area just south of Shellharbour through to Bega (~157km north of Victorian boarder) [35].After the fires, the south-east corner of NSW experienced some of the largest percentage increases in erosion due to reduced soil stability-a common consequence of bushfires [35].Bushfires can rapidly impact estuaries through such erosion, especially when coupled with post-fire rainfall, and can thus influence estuary hydrology, turbidity and nutrient inputs via sedimentation.Benthic habitats including seagrass beds and shellfish reefs [28].whichare important are important habitats for bream [76,77], can be affected by this sedimentation [28].The intensity of the bushfires in the south-east may have therefore reduced the habitat available to bream and contributed to their mortality in this region to a greater extent, possibly resulting in the lower catches during 2019/20.Silva, Doyle [78] identified A. butcheri (black bream) during assessments of aquatic biota mortalities within South Coast waterways following the 2019/20 fires.It is also possible that the change in water quality simply led to species such as bream avoiding fire-impacted estuaries in the South Coast, thus leading to lower catches.For ocean caught bream, total catch in the Central Coast was lowest 2013/14 with moderate to high recreational catches estimated for 2017/18 and 2019/20, respectively (Table 7; point ii).In the South Coast however, the oceanic catch of bream was similarly high in 2013/14 and 2017/18 and lowest in 2019/20 (Table 7; point ii).The drought conditions recorded between 2018 to early 2019 resulted in reduced water in the soil, rivers and groundwater across much of Australia but had contrasting effects on streamflow from south to north in eastern Australia [79].These conditions are one of many possible factors that may have contributed to the differences observed in oceanic bream catch.In NSW, the drought was shown to decrease coastal streamflow more severely in the Central and North Coast compared to the South Coast and may partly explain the opposing patterns observed for oceanic bream catch in the Central and South Coast [79].This north-ward decline in streamflow is likely to have increased estuarine salinity in the Central and North Coast and could have triggered downstream migration into oceanic habitats for marine-estuarine opportunist species like bream [29,80].Accordingly, Gillson, Suthers and Scandol [29] showed that NSW commercial landings of A. australis (yellowfin bream) were greater in oceanic habitats during drought due to an increase in their abundance and thus catchability following estuarine emigration.A similar phenomenon may have occurred with oceanic recreational landings in the Central Coast following the drought period of 2018 to early 2019.

Total
Although the interaction between region and year was not deemed important in explaining variation in fishing effort in either estuarine or oceanic waters, independent bioregional and temporal effects were found.When temporal effects were examined without accounting for bioregional effects, 2019/20 consistently showed increased levels of effort when compared to 2017/18 (Table 7; point iii).Based on trends reported for the state's freshwater recreational fishery and in other parts of Australia [39,50], it was expected that adverse fishing conditions caused by the bushfires and flooding, in particular, might cause a decline in coastal effort by 2019/20.In contrast, the high effort observed for 2019/20 here may reflect the early influence of socio-economic changes caused by the COVID-19 pandemic on coastal recreational fishing, over and above the impacts of the extreme weather events.As governments tried to control the spread of COVID-19, measures were taken in early 2020 (and enforced from March in Australia) that restricted human movement, led to temporary cessation of non-essential activities, such as tourism, and slowed many local economies [81,82].In NSW, these changes resulted in high proportions of coastal metropolitan populations working from home and other portions experiencing job losses or reduced work hours.With less time committed to work or commuting, fishers in NSW may have increased their engagement with the relatively safe (i.e.socially isolated) and low-cost outdoor activity of recreational fishing, as was reported in other parts of the world where fishing activity exceeded pre-pandemic levels [38,83,84].
For estuarine and oceanic snapper catch as well as overall CPUE, inter-annual variation was not detected, providing little evidence that these metrics were influenced by the extreme weather events for this species.Adult snapper are generally found in offshore oceanic waters [85], where environmental perturbations caused by events such as bushfires, droughts and floods are likely to be of less consequence.Although extreme weather events can impact the estuarine and near-shore nursery environments inhabited by juvenile snapper [85], past research suggests that juvenile snapper are well-adapted to the shifts in salinity, turbidity and oxygen availability [86][87][88] brought about by extreme events.Snapper may have therefore had a greater resilience and likelihood of recruiting into recreational fisheries in both estuarine and oceanic environments despite the extreme weather events.CPUE could only be confidently estimated across all species caught using line-based methods (~92 individual species/ taxa).The inclusion of such a large number of species with divergent patterns of habitat-use, biology and differing responses to factors like extreme weather events may have introduced substantial noise into the catch component of CPUE and therefore made it challenging to detect interannual variation.Overcoming this limitation in the future could involve quantifying species-specific CPUE through customized surveys that do not rely on broad assumptions about targeting behaviour.

Conclusion
This study shows interesting variation in fishery metrics in a coastal recreational fishery that may be driven by socio-ecological dynamics.It also highlights the utility of off-site surveys of recreational fishers in monitoring this variability.As climate change and related extreme weather events continue to affect the world's coastal environments, fisheries are likely to be impacted and will require ongoing adaptive management to ensure their sustainability.Quantitative monitoring studies can contribute critical data to assist in identifying problematic trends and in making projections.This, in turn, facilitates the development of targeted management strategies to address future threats to recreational fisheries [89,90].Such strategies may include habitat enhancement through the deployment of purpose-built artificial reef systems in estuaries affected by sedimentation to protect bream populations in regions prone to bushfires and floods [91].Another targeted management strategy could involve the introduction of water management plans that determine and implement environmental flow requirements specific to estuarine ecosystems and manage these flows holistically from catchment to coast [92].This strategy might ensure adequate stream flow for fish in estuaries within drought affected regions.

Fig 1 .
Fig 1. Map of New South Wales (NSW), Australia indicating the relative sizes of the North Coast, Central Coast and South Coast regions.The state boundary shown in the map was generated by Geoscience Australia as part of the Australian Marine Spatial Information System (https://amsis-geoscience-au.hub.arcgis.com/)and the Bioregion boundaries are generated from the Integrated Marine and Coastal Regionalisation of Australia v4.0 (https://fed.dcceew.gov.au/datasets/erin::integrated-marine-and-coastal-regionalisation-of-australia-imcra-v4-0-provincial-bioregions/about).The source maps are released under the Creative Commons Attribution (4.0 and .0respectively) International Licence.Included is a timeline showing the sequence of extreme weather events that took place within NSW during the 2017/18 and 2019/20 recreational fishing surveys.Note that a preceding 2013/14 survey is not included in the timeline.https://doi.org/10.1371/journal.pone.0305106.g001

Table 6 .
Results of the zero-inflated Poisson generalized linear mixed effects models used to examine the influence of year (3 levels: 2013/14, 2017/18 and 2019/20) and region (3 levels: North Coast, Central Coast and South Coast) on the catch per unit effort (CPUE; numbers/fisher days) estimated for line-based recreational fishing activity in the estuarine and oceanic waters of NSW during 2013/14, 2017/18 and 2019/20.

Table 3 . Results of the generalized linear mixed effects models used to examine the influence of year (3 levels: 2013/14, 2017/18 and 2019/20) and region (3 levels: North Coast, Central Coast and South Coast) on the estimated total number of bream and snapper caught (retained and released) through line-based recreational fishing activity in the estuarine and oceanic waters of NSW during 2013/14, 2017/18 and 2019/20. The
models, which assumed a negative binomial distribution, applied the Wald test (α = 0.01) to examine the null hypothesis that a β i = 0, where β i s were the regression coefficients for region, year or their interaction.The total sample sizes (n = number of long-term Recreational Fishing Licence holders) of each model are also shown.https://doi.org/10.1371/journal.pone.0305106.t003

Table 4 . The results of Tukey's post-hoc tests (α = 0.01) comparing different levels of year (3 levels: 2013/14, 2017/18 and 2019/20) or region (3 levels: North Coast, Central Coast and South Coast) in terms of the total catch of bream, total catch of snapper, total effort and total catch per unit effort (CPUE) as estimated for line- based recreational fishing in the estuarine and oceanic waters of NSW during 2013/14, 2017/18 and 2019/20.
https://doi.org/10.1371/journal.pone.0305106.t004was significantly greater in the Central Coast compared to both the South Coast and North Coast (Table 4, Fig 3A).Temporally, the only significant difference in effort was between 2019/20 and 2017/18 with the former year having the greatest effort (Table 4, Fig 3A).

Table 7 . Summary of the main inter-annual differences in fishery metrics observed in this study and a brief description of the socio-ecological factors that may explain each observation. Point Inter-annual difference Potential explanatory socio-ecological factor
https://doi.org/10.1371/journal.pone.0305106.t007